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This article explores the measurement of temperature in transient states, 
utilizing the principles of heat transfer and thermal-electrical metaphor. The 
study focuses on the nonlinear thermal resistances present in various locations 
within a distribution transformer, while taking into account variations in oil 
physical variables and temperature loss. Real-time data obtained from heat 


run tests on a 250-MVA-ONAF cooled unit, conducted by the transformer 
manufacturer, is used to verify the thermal designs. The observations are then 
compared to the loading framework of the IEC 60076-7:2005 system. The 
findings of this research provide a better understanding of temperature 
measurement in transient states, particularly in distribution transformers, and 
can be applied to the design and development of more efficient and reliable 
transformer systems. 
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1. INTRODUCTION 

The significant proportion of infrastructure investment in generation and distribution power stations 
is characterize by power transformers. Moreover, the transformer is one of the costliest elements of an energy 
grid as the power transformer disruptions have a significant economic effect on the operation of the 
transmission grid. Identifying their status is thus important for achieving the objectives of optimizing return on 
investment and reducing the overall costs related to the operation of the transformer. The warm heating rate is 
amongst the most model aspects controlling the average lifespan of a transformer. Among the most important 
metrics in evaluating the life of power transformer is the winding warm temperature, as the highest aging rate 
happens at the highest point that encounters the greatest temperature. In addition, for the transformer to have a 
healthy average lifespan, the warm temperature must blow up the allowed maximum amount. 

It is commonly acknowledge that for a 6 °C rise in temperature, the insulation loss roughly doubles 
[1]. An effective warm temperature estimation is thus critical for both suppliers and users. There are a few 
ways to calculate the temperature; the first is to use fiber optical heating elements located at the windings' that 
anticipated the warmth. Thermal sensors mounted to the bottom of the fiber optics are normally mount between 
the isolated connector and the sprocket while their impulses are emit from the tank through the optical fiber. 
From over decades, major changes resolve the fiber optic being far too brittle and sensitive handling is need. 
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Precise calculation of warm temperature can be achieve from using this process. The IEEE and IEC loading 
manuals [2], [3] used to measure the temperature of the warm using heat-run experiment data and analytical 
variables. In addition, the temperature of the transformer rely on the exact failures in the transformer produced. 
This implies that the different transformer root and wind losses are important for deciding the winding intensity. 
The IEEE C.57.120 [4], [5] models can measure transformer losses. Even so, the real losses are required to 
determine the actual wind level, for example, the losses should be measure based on data. Iqteit and Yahya [6] 
described an automated transformer loss calculation, which requires a comprehensive measurement operation. 
Fuchs et al. [7] also constructed an enhanced real time tracking system for transformer damages, differentiating 
among iron core and copper losses. In addition to measuring the iron core, researchers incorporate eddy current 
and hysteresis loses and then use ohmic wind and stray losses to assess the conduction loss. Nevertheless, 
another such distinction among ohmic wind and stray losses is important since the stray losses arise in various 
parts of the transformer [8] and hence have only a slight influence on the intensity of the winding. 
Consequently, the temperature is primarily influence by losses in ohmic winding [9]-[11]. 

Figure | provides a simple thermal management system for power transformers, for which oil intensity 
within the windings is supposed to rise linearly from start to end, while the temperature differential here 
between pipe oil and the winding coil is steady during the winding phase [12], [13]. The change in the warmth 
temperature is greater than the change in the concentration of the coil at the tip of the wind, which can see in 
Figure | whenever the changes due to stray losses are factor into the equation. Even so, this conventional 
method of measuring temperature s has found to be insufficient in accordance with the real data recorded by 
fiber optic sensors [14], [15]. Blume et al. [16], [17] noted that there's still a time difference between the rising 
in top-oil temp as well as the rising in duct oil temp under dynamic loading. The consequence of this effect is 
that the wind hotspot rate is greater than the “IEEE loading guide clause 7” method, which is expected. 

It was introduced in [16] that calculations to accommodate for the oil duct temps, the difference in 
stator winding with temp, the difference in oil viscosity, the impact of the tap direction, as well as the difference 
in atmospheric temp (temp) as during load period, that were overlooked in the system of “Clause 7”. These 
updated calculations have been used in the “IEEE loading guide” [3] as option to the hotspot temperature 
measurement procedure termed as the “Annex G method”. In this article, a quantitative analysis was perform 
to measure the hotspot temperature amongst the two techniques specified in the IEEE loading guide [3]. In key 
findings of continuous load, load demand pattern and short-period overload scenarios, four transformers with 
various cooling mechanisms of “oil natural air natural (ONAN), oil natural air forced (ONAF), oil forced air 
forced (OFAF), and oil directed air forced (ODAF)”, the measurements of whom are provide in the appendix 
section, are being used to test the two approaches. Events of ONAF and OFAF measurement have been 
confirmed by thermal experiments and [16], [17] calculations. A MATLAB/simulink software was used to 
analyze the hotspot temperature and the consequential loss-of-life. 
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Figure 1. Basic thermal diagram inside the transformer tank [6] 


2. TRANSFORMER THERMAL MODELS 

The sources from [4]—[8], and [18]-[24] discuss earlier stuff performed in this field of transformer 
thermal modeling. For all the temperature measurements and wind hotspot forecasts, which have been the IEC 
and IEEE scientific formulae, the normal approaches followed [18], [19] are inadequate for today's industry 
competency [20], [21]. Furthermore, basic temperature rise predictions, which have been conducted even now 
days by the traditional scientific formulae, are not just useful indicators for direct temp calculation or emulation. 
In an approach to strengthen the precision, further advanced models are now being updated with IEEE and IEC 
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loading guidelines aimed at a fair analysis of the oil temp within the winding, taking into account changes in 
winding resistance, oil viscosity and oil inertia. The energy balance formula will describe a thermal process. 


qx dt = Cy x do + ame x dt (1) 


The given equation includes several variables such as g, Cn, 8, Rth, and Oamp. These variables represent 
the heat generation, thermal capacitance, temperature, thermal resistance, and ambient temperature, 
respectively. The formula can be interpreted as (2). 


= a0, 9-Pamp 
q=CihX 7+ a (2) 
Then, when the describtion a small electronic RC circuit as shown in Figure 2 makes a formula based on both 
the first principle of Kirchhoff and Ohm. 

(3) 


: du u 

l= Coy x at + Rel 

The (3) consists of four distinct variables, i, C,;, Re,, and u, which represent key components of the 
electrical system, including electrical current, electrical capacitance, electrical resistance, and electrical 
voltage, respectively. Merely, we achieve the comparison among electrical and thermal mechanisms by 
contrasting (3) and (2) as shown in Table 1. Figure 3, the analog thermal circuitry for the electrical process has 
provided. 


Lt fd | 


Figure 2. An electrical RC circuit Figure 3. The analogous thermal circuit 


Table 1. Thermal-electrical analogy 


Thermal Electrical 
Generated heat q Current i 
Temperature 6 Voltage u 
Resistance Rin ‘Resistance —- Re, 
Capacitance Cyn Capacitance Cy, 


2.1. The non-linear thermal resistance 
The nonlinear oil thermal resistance, Ryp—oi, (m?K)/W, in accordance with the findings of energy 
transfer the corresponding formula is express in (4). 


— 1 _ Abit 
Rio > ang (4) 


In the (4), / represents the heat transfer coefficient, A is the area, A@,;, corresponds to the oil temperature 
gradient, and g indicates the amount of heat generated by the relevant losses. The natural flow of convection 
oil across longitudinal, leaned and vertical sheets and tubes can be determined by the following integral 
approach on the basis of heat transfer theory [14]-[16]. 

N, = C x [G, x P,]” (5) 
The coefficients C and n depend on whether the oil circulation is laminar or turbulent, and they are obtained 


empirically. The Nusselt number (N,,), Prandtle number (P.), and Grashof number (G,.) describe in (3) to (6). 
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me 6) 
Coilx 

B= (7) 
13 xp? .XgXBX (AO it) 

G, = ExPiuroxPs(aeon (8) 


Where the L represents the characteristic dimension of the system, which could be the length, width, or 
diameter of the transformer. The symbol g stands for the gravitational constant, k denotes the thermal 
conductivity of the transformer oil, po; is the density of the oil, 6 represents the oil thermal expansion 
coefficient, C,;; 1s the specific heat of the oil, 1 represents the viscosity of the oil, and A@,;) stands for the oil 
temperature gradient (K). The continuity equation is achieve by supplementing the (6) to (8) into (5). 


Wet Cx | (Saat) x (ee 0) 
k k pe? 

Including all transformer isolation oils, it's indeed usually correct that the temp difference of the 
operating temp is far greater than many of the other oil variables [18]—[20], [25]. Itis also possible to substitute 
all physicochemical characteristics of the oil except perhaps the viscosity in (9) with a standard. Even so, if the 
effect of all oil variables needs to be weight, the foregoing actions should take. The (10) for the thermal 
conductivity will be estimate by the following: 


: (a=n)\ ™ 
h=C,x (20 x atte) (10) 


where C1 is assumed to be a constant, and is now expressed as (11). 


3n-1\,() 
C= 0x [g xi] (11) 
As well as the temp difference of all oil variables is shown by the corresponding equation [20], [21]. 
lait 
w=A,Xe Ooi +273 (12) 
Coit = Az X Ag Ooi (13) 
Poit = As X Ae Aoi (14) 
k = Az X Ag9oi1 (15) 


B= Ay (16) 


Table 2 lists the nine conditions for the two structural component oils. 


Table 2. Insulation oil constants [18], [20] 


Oil/constant — Transformer oil Oil/constant Transformer oil 
Al 0.13573x10° A6 -0.659 
A2 2797.3 AT 0.124 
A3 1960 A8& -1.525x10-4 
A4 4.005 A9 8.6x10-4 
AS 887 


2.2. The top-oil temp model 
As a thermal chain, Figure 4 offers the top-oil temp model. Centered on the principle of heat flow and 
thermal-electrical analogy [1]-[3], [14]. Where, the total heat generated by losses in a transformer is 
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represented by q;o¢, with qe indicating the heat generated by no-load losses and q, representing the heat 
generated by load losses. The equivalent thermal capacitance of the transformer oil is denoted by c¢j_9;; , and 
Ooi, tepresents the top oil temperature. The non-linear oil to air thermal resistance is represented by 
Ren—oit—air, 20d Ogmp represents the ambient temperature. Two basic heating elements reflect the heat 
produced by certain no-load and load transformer falls and the atmospheric temp is described as a value 
obtained of temp [1], [2]. The nonlinear system seen in Figure 3 for the thermal system is: 


doit , (Goit—Pamb) 


+ q; = Crn—oi X + 17 
Ife T 1 th-oil * “Ge Raccasae ce) 
if we replace the non-linear heat transfer formula (4) with (17), the corresponding model is solved. 
db oi (8oi1-9 ) 
dge + 1 = Chon X St (18) 


hxA 


After which, by replacing the thermal conductivity variables for formula (10), h, the differential equation is 
modified to, 


n 


n 
H u 
p2xBXxCoipxkKA-M/N} —_\ p2xBxcgipxkG-m/n 


(We +) X ea “yo ee Crnoit X weit + (96:1 — Pamp)**"(19) 
after that, the parameter can be described as (20)-(24). 

LH = pu X Uratea (20) 

P = Ppu X Prated (21) 

B = Bou X Bratea (22) 

C = Cpu X Crated (23) 

k = kyu X Kratea (24) 


As well as the parameters that follow, the measured thermal non-linear instance, R¢p_oi1—airratea, the 
measured temp of the top oil rises above the ambient temp, A@ i -atea: 


AO oitratea = (Fe + a) Reh-oil-airrated (25) 


x 
rated 


The rated top — oil time constant, t_(oil, rated). 
Toil,rated = Rth-oil-air,rated, x Cth-oilrated (26) 
The proportion of load losses to no-load losses at nominal current, R, 


yacde (27) 
fe 


and the load factor, K. 


I 


= 28 
Trated ( ) 
Then expression (19) is simplified to its exact from (29). 
4 n n 
1+R+K Hpu Hpu Foi 
x = xX Aa; = = X Toi x + 
1+R (eae wm) oil,rated (=e c=] oilrated dt 
(9oi1-Pamp)**” 29 
net (29) 


oilrated 
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2.3. The hot-spot temp model 

The comparison of non-linear opposition here between protective coating surface and the wind 
surface, the hotspot temp method, which would be centered on the top oil temp will be established. The 
framework is founded on the principle of traditional energy transfer. The hot-spot temp model is also depicted 
as a spectral network, comparable to the traditional heat transfer practice for the top-oil temp framework and 
the non-linear heat capacity as Figure 5 shows [1], [2]. 


Kth-oil- 
tot Gait theater hs Reiehs-oil 


Ve qi a le Gdn | Ctn-wan ie 
p) qd | p si if 


Figure 4. The top-oil temp model Figure 5. The hot-spot temperature model 


+ 
+ 


The parameters in the figure include qwq@y, which represents the heat generated by winding losses, 
Ctn—wan, Which is the thermal capacitance of the winding, @,,, which is the hot-spot temperature, Rrp—ns—oit, 
which is the nonlinear thermal resistance between the winding and the oil, and 6,,;,, which is the temperature 
of the top oil. The vibrational wind is determined using an expression for oil thermal conductivity. 


Ren—ns—oit = Rten-wan + Ren—inswt + Ren—insut—oit (30) 


The thermal resistances used in the model include the winding thermal resistance Ryp_wan, Winding insulation 
thermal resistance R¢p_insy;, and non-linear winding insulation to oil thermal resistance Ryp_insui—oi. BY 
comparing the resistances provided in (30), the properties of the transformer can be analyzed: 


Ren-insut-oit >> Rten-wan (31) 
Ren-insut-oit > Ren—insut (32) 


for the warm conditions recorded on the insulation's outermost layer twisted across the conductors [22], [23]. 
The critical value for the nonlinear wind to thermal conductivity of oil is therefore: 


1 
Ren-ns-oil = Faq (33) 


XA 
for both the top oil temp framework, the expression (27) is identical to (4), so the coefficient for the temp 
profile, h, is completely equivalent to the temp profile in (10). 


3 (a-n)\" 
h= C, x (20 x Hate) (34) 


In which all the parameters of the oil are reassessed at the top temp of the oil and A@,,, is now the hot- 
spot to top-oil temp gradient. The nonlinear system is shown Figure 4 for the thermal circuitry is: 


dOns  (Ons—9oi1) 
dwan = Cen-wan X +o (35) 
dt = Rth—hs—oil 


if the non-linear thermal conductivity model replaced (34), with (35), the corresponding model is solved. 


dons 


dwan = Cen-wan X ts + Caso (36) 


hxA 


After which, by replacing the thermal conductivity expression, (34) with (36), the dynamic model is 
modified to: 
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dbps 
at 


n n 
fa K 
(ae ea) = mae 


a yl+n 
C1xA oT + (Bps — Ooi)? *"(37) 


Iwan * X Crh-wan * 


The average hot-spot non-linear to top-oil thermal conductivity, Rtn—ns—oil,ratea» The measured hot-spot 
temp rises above the temp of the top oil, A@n57atea, aS Shown in (38). 


A®hns rated, = Awdn rated x Rth-hs—oil,rated =H x Gr (38) 


Whereas the hot-spot factor H and the rated average winding to average oil temp gradient g, are defined in 
reference [24]. The rated winding time constant Tyanratea, 1S: 


Twdnrated, as Rth—hs—oil,rated* Crn—-wanrated (39) 


Pnst Ox “sated Ph) 


+ x ( 40 
wot) Peddy,pu Ost OK ( ) 


Dijana eis) = Pdcpu * ( 


where, Pac pu(Pns) and Peaaypu(Pns) The activity of the DC and eddy drops is defined as a temp function, the 
DC losses vary with the temperature directly, while the eddy losses vary with the temp inverse proportion. 0, 
It is the loss adjustment temperature variable, equivalent to 225 for aluminum and 235 for copper. It implies 
that the definitive formula is: 


n n 
Upu Upu d8ns _, (Ons—9oit)*” 
k? X Pwanpu(Ons) } x | =———*"* aa) APhs rated = = X Ths,ratea X opanaitaet 
{ sath tes . j PB uXBpuXCoilpuXks,, Bie aaa PB uxBpuXCoilpuxk yy, mie Se ore dt APhsrated 


the effect of capacitors that may optimize the maximum electrical power saving is shown in [26], where [27] 
calculates the transformer minimum energy loss. 


3. SIMULATION MODEL 

In MATLAB/simulink, the thermal coefficients are patterned. The coefficients are solved 
simultaneously using the system of Runge-Kutta. The requisite variables for the model are based on data 
obtained from the standard heat run test carried out by the supplier of the transformer. Each one phase the 
formula of the top oil model (29) is resolved by K, Upy, Ppw Bow Coit pw Kpuand Oamp are the input variable 
while 6,;; The output vector, although both are time constants, is t. Figure 6 presents a schematic diagram of 
the thermal model for measuring temperature during transient states in a distribution transformer. The applied 
top oil design is demonstrated in Figure 6(a). 

As seen in the condensed Figure 6(b), the measured top oil temp is the input atmospheric temp for the 
hot spot method. Each one phase the formula of the top oil model (29) is resolved by K, 
Hpw Pow Bow Coipw Kpw Pwanpu(hs) and 6,;; are the input variable while @, the output vector, although both 
are time constants, is t. And use the block diagram seen in Figure 7, the applied top oil design is demonstrated.is 
the output variable which all are functions of time, t. 


n 
1+R+K2 Upu 
Where Q1 = Tear x ——_— > en) x AO (it rated) 
PouXBpuXCoitpuxkpy” 
1 
Q2 — n 
fe (i—n) XT (oil,rated) 
pe uXBpuXCoilpuxkpu 
n 
m 
=e ge2 Pus = 4s =. 

Q3 = k* x Pwanpu(9ns) x aay} * AG (hs,rated) 


2 
Pou*P puX CoitpuX Kou" 
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Hot-spot ; 
model 


(b) 


Figure 6. The schematic diagram model: (a) block diagram of top oil model and (b) simplified diagram of 
the thermal model 


Figure 7. Block diagram of hot spot model 


4. SIMULATION RESULATS 

This section presents the temperature results obtained from three different transformer modules and 
various types of tanks during various load tests. These results were obtained using equations (29) and (41) and 
were compared to the IEC 60076-7:2005 loading guide method. The data used in this analysis was derived 
from the normal heat training session conducted by the transformer company. 


4.1. Transformers with external cooling 

Winding temperature restricts the loading of the transformer, so the temperature of the power 
transformer must maintain beyond a certain boundaries prescribed by regulatory requirements for maximum 
load and normal ambient temperature. The inrush current temperature is not standardized, and the decentralized 
control is usually the winding's elevated temperatures portion, named the winding hot spot temperature. The 
insulated temperature is the key element in the ageing of the transformer. The sheet separation is translocated 
with the temperature and time reflecting the end of the insulated materials' operation, which is known as the 
terminal phase of the transformer. 


Table 3. The load steps for the 250 MVA transformer 


Time period (minutes) _ Load 


0.0-187.4 1.0 
187.4-364.9 0.6 
364.9-503.4 1.5 
503.4-710.0 0.3 
710.0-735.0 241 
735.0-750.0 0.0 


Figure 8 shows the rated impedances for the 250 MVA transformer, which were 23048 x 1.5 
percent/118/21 kV. Figure 8(a) and Figure 8(b) illustrate the load steps for the transformer as presented in 
Table 3. The oil flow through the 118 kV and 230 kV windings was directed in a zigzag pattern by oil directing 
circles. The top-oil temperature is presented in Figure 9. 
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Figure 8. Load steps and hot-spot temperature of the 250 MVA ONAF-cooled transformer: 
(a) load steps and (b) hot-spot temperature of the 118 kV winding 
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Figure 9. The top-oil temperature of the 250 MVA ONAF-cooled transformer 


5. CONCLUSION 

This article analyzes an effort by the established thermal electric comparison approach to analyze the 
specific hot-spot temperature physical analysis for more detailed temperature measurements during load 
variations. The specific thermal model for the top-oil is often analyzed by the thermal electrical analogy 
approach used, taking into account the influence of atmospheric temperature on the top-oil temperature. The 
constants (n) used to define the spectral gradient for the top oil numerical simulation model and the hot spot 
thermal model are crucial for different cooling mechanisms and transformer architectures. These models were 
defined in the equations used throughout this research. In the numerical simulation, both oil variable 
modification and loss variance with temperature are taken into consideration. The findings acquired by the 
numerical simulation take into consideration of all variations in oil physical variables and the difference in loss 
with temperature is in fair accordance with the findings plotted by the thermal model. The results stated that, 
it is worth noting that the IEC 60076-7:2005 loading guide method results less well with the results obtained 
by thermal model. 
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